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Abstract 

We reinvestigate how generic off-diagonal cosmological solutions depending, in general, on 
all spacetime coordinates can be constructed in massive and f-modified gravity using the anholo- 
nomic frame deformation method. There are constructed new classes of locally anisotropic and 
(in) homogeneous cosmological metrics with open and closed spatial geometries. By resorting 
such solutions, we show that they describe the late time acceleration due to effective cosmolog¬ 
ical terms induced by nonlinear off-diagonal interactions, possible modifications of the gravita¬ 
tional action and graviton mass. The cosmological metrics and related Stiickelberg fields are 
constructed in explicit form up to nonholonomic frame transforms of the Friedmann-Lamaitre- 
Robertson-Walker (FLRW) coordinates. The solutions include matter, graviton mass and other 
effective sources modelling nonlinear gravitational and matter fields interactions with polariza¬ 
tion of physical constants and deformations of metrics, which may explain dark energy and 
dark matter effects. However, we argue that it is not obligatory always to modify gravity if 
we consider effective generalized Einstein equations with nontrivial vacuum and/or non-minimal 
coupling with matter. Indeed, we state certain conditions when such configurations mimic in¬ 
teresting solutions in general relativity and modifications, for instance, when we can extract 
the general Painleve-Gullstrand and FLRW metrics. In a more general context, we elaborate 
on a reconstruction procedure for off-diagonal cosmological solutions which describe cyclic and 
ekpyrotic universes. Finally, there are discussed open issues and further perspectives. 
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1 Introduction 

The "ekpyrotic" and "new ekpyrotic" mechanisms and cyclic models have been elaborated as 
alternatives to standard big bank inflationary cosmology |2]. Another alternative came from the 
idea that graviton may have a nontrivial mass as it was proposed in the Fierz and Pauli work 
[I]. For recent reviews and related f(R) modifications and cosmological models, in general, with 
non-minimal coupling and dilaton-brane cosmology, local anisotropies and/or effective modelling of 
massive gravity, see respectively p.m and [5]. The modern version of a ghost free (bimetric) massive 
gravity theory were made in a series of papers (recent reviews can be found in |6]); Developing generic 
nolinear versions of the Fierz-Pauli theory, the so-called vDVZ discontinuity problem was solved via 
the Vainshtein mechanism [7j (avoiding discontinuity by going beyond the linear theory). There are 
also more recent approaches based on DGP model |8]. During a long time, none solution was found 
for another problem with ghosts because at nonlinear order in massive gravity appears a sixth scalar 
degree of freedom as a ghost. This problem was cosndiered in a paper by Boulware and Deser [9] 
together with similar issues related to the effective field theory approach. 

A considerable amount of work has been made in order to undestand the implications and hnd 
possible applications of of such ghost-free models. The most substantial progress was made when de 
Rham and co-authors had shown how to eliminate the scalar mode and Hassan and Rosen provided 
a complete proof for a class of bigravity / bimetric gravity theories, see nm. In such approaches, the 
second metric describes an effective exotic matter related to massive gravitons and does not suffer 
from ghost instability to all orders in a perturbation theory and away from the decoupling limit. 

It is important to mention that the first black hole solutions in a nonlinear massive theory 
were found in the context of high-energy physics and for off-diagonal and/or higher dimensions 
generalizations m- A general decoupling of generalized Einstein equations can be proven following 
the anholonomic frame deformation methods, AFDM, jl2]. Here we also note that the possibility 
that the graviton has a nonzero mass fl results not only in fundamental theoretical implications 
but give rise to straightforward phenomenological and cosmological consequences. For instance, a 
gravitational potential of Yukawa form ~ e“^^/r results in decay of gravitational interactions at 
scales r > /i“^ and such an effect may result in the accelerated expansion of the Universe. This way, 
a massive gravity theory, MGT, provides alternatives to dark energy and, via effective polarizations 
of fundamental physical constants may explain certain dark matter effects. We can treat this as a 
result of generic off-diagonal nonlinear interactions. Recently, various cosmological models derived 
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for ghost free (modified) massive gravity and bigravity theories have been elaborated and studied 
intensively |T?11?1 IH]. 

The goal of this work is to construct generic off-diagonal cosmological solutions in MGT and state 
the conditions when such conhgurations are modelled equivalently in general relativity (GR). It ex¬ 
tends the letter variant of Ha with new results and more details on proofs of results and constructing 
exact solutions and further analysis of results and speculation on reconstruction mechanism. 

As a hrst step, we consider off-diagonal deformations of a "prime" cosmological solution taken in 
general Painleve-Gullstrand (PG) form, when the Friedman-Lamaitre-Robertson-Walker (FLRW) 
can be recast for well-dehned geometric conditions. Such constructions are performed in section [2j 
At the second step, the "target" metrics will be generated to possess one Killing symmetry 
(or other none Killing symmetries) and depend on timelike and certain (all) spacelike coordinates. 
In general, such off-diagonal solutions are with local anisotropy and inhomogeneities for effective 
cosmological constants and polarizations of other physical constants and coefficients of cosmological 
metrics which can be modelled both in MGT and GR. We consider the method of constructing 
generic off-diagonal cosmological solutions in section |3l 

Then (the third step), we shall emphasize and speculate on importance of off-diagonal nonlinear 
gravitational interactions for elaborating cosmological scenarios when dark matter and dark energy 
effects can be explained by anisotropic polarizations of vacuum and/or de Sitter like conhgurations. 
We provide examples of off-diagonal solutions with solitonic conhgurations, see section 01 

A reconstruction mechanism for oh-diagonal cosmological solutions with modihed gravity and/or 
massive graviton ehects is elaborated in section [51 Finally, we conclude the paper in section [6l 

2 Equivalent modelling of /—modified and massive gravity the¬ 
ories 

We shall work with MGTs modelled on a pseudo-Riemannian spacetime V with physical metric 
g = and hducial metrics. On massive gravity, see reviews (Hj and on geometric methods in 

gravity and constructing exact solutions, see Refs. na . In addition to well-known approaches with 
diadic and tetradic (vierbein) variables, we shall work with nonholonomic manifolds when certain 
classes of frame transforms can be adapted to preserve a chosen splitting of the nonlinear and linear 
connection structures into some standard components (for instance, dehning the Levi Givita, LG, 
connection) and distortion tensors which can be hxed to be zero if additional constraints are imposed. 

Using diherent nonholonomic frame variables, the action for our model can be written in two 
forms. 


S = K„;,) + "L] (1) 

= ”ri- (2) 

In above formulas, the physical/geometrical objects with "hat" and/or written in "boldface" form 
are considered for a conventional 2-1-2 splitting when the two dimensional horizontal, h, and the two 
dimensional vertical, v, coordinates are labelled (respectively) in the form = (x®, ?/“), or u = {x,y), 
with indices i,j, k... = 1,2 and a, b ,... = 3,4. The scalar curvature R is for the LG-connection V. 

We write by R the scalar curvature for an auxiliary (canonical) connection D uniquely determined 
by two conditions: 

1. it is metric compatible, Dg = 0, and 


3 






2. the h- and n-torsions are zero (but there are nonzero h — v components of torsion T completely 
determined by g) for a conventional splitting 

N ; TV =hV®vY, 

with local coefficients of the so-called nonlinear connection, N-connection structure, labeled in 
the form N = {iVf}. 

General frame transforms can be parameterized in the form = A°‘'^{u)da>■, where the matix 
is non-degenerate in a hnite, or inhnite region of V and da' = djdvA . Using such we can always 
re-de£ne the geometric and physical object with respect to a class of N-adapted (dual) bases 

ea = (ej = di- N-db, Ca = da = d/dy'") and 

= dx\ = dy^ + N^dy^), (3) 

which are nonholonomic (equivalently, anholonomic) because, in general, there are satished relations 
of type 

for certain nontrivial anholonomy coefficients fU2^(M). The Einstein summation rule on repeating 
indices will be applied if the contrary is not stated. 

The connection D allows us to decouple the held equations in various gravity theories and con¬ 
struct exact solutions in very general forms. Here we note that the distortion relation from the LC 
connection 

D = V + Z[f] 

is uniquely determined by a distorting tensor Z completely dehned by T and (as a consequence for 
such models) by (g, N). The main idea of the AFDM [12] is to use D as an auxiliary [for the (pseudo) 
Riemennian spacetimes] and/or [with nonholonomically induced torsion] connection which positively 
allows us to decouple gravitational held equations for very general conditions with respect to N- 
adapted frames (|2]). We can not perform such a decoupling in general form if we work from the very 
beginning with the data (V,c?q') but there are proofs that this is possible for (D,eQ;g,N). Having 
constructed integral varieties (for instance, locally anisotropic and/or inhomogeneous cosmological 
ones), we can impose additional nonholonomic (non-integrable constraints) when —)■ V and 

R ^ R, where R is the scalar curvature of V, and it is possible to extract exact solutions in GR|/| 
The MGTs with actions of type ([ID generalize the so-called modihed f{R) gravity and the ghost- 
free massive gravity [I^. We shall follow some conventions from m. There will be used the units 
when h = c = 1 and the Planck mass Mpi is dehned via Mpi = l/SvrG, with 4-d Newton constant 
G. We write 5u^ instead of d'^u because there are used N-elongated diherentials as in (|2D- A model 
can be specihed by corresponding constants, dynamical physical equations and their solutions in 
corresponding variables. It will be considered that y = const is the mass of graviton. For LG- 
conhgurations, we can hx conditions of type 

fiR) - ^W(g^„ Kap) = /(R), or fiR) = /(R), or /(R) = R, (4) 

which depend on the type of models we elaborate and what classes of solutions we wont to construct. 
Such conditions can be very general ones for arbitrary frame transforms and N-connection deforma¬ 
tions. We emphasize that it is possible to hnd solutions in explicit form if we chose the coefficients 

^There will be considered also left and right up/low indices as labels for some geometric/physical objects. 
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{A^“} and the local frames for D when R = const in such forms that daf{R) = (d^f) x daR = 0 
but, in general, daf{R) 7 ^ 0 . 

The equations of motion for our nonholonomically modihed massive gravity theory can be written 

- ^fiR)g,u + = M-f(5) 

where Mpi is the Plank mass, is the Einstein tensor for a pseudo-Riemannian metric and 
D, is the standard matter energy-momentum tensor. We note that in above formulas, for 
D —)■ V, we get —)■ with a standard Ricci tensor R^^ for V. Such limits give rise to original 
modihed massive gravity models but it is important to hnd some generalized solutons and only at 
the end to consider additional assumptions on limits and nonholonomic constraints. It should be 
emphasized that the effective energy-momentum tensor in ([5]) is dehned by the potential of 
graviton U = U 2 + a.-M'i + ^ 4 ^ 4 , where 03 and Q !4 are free parameters. The values U 2 MZ and W 4 are 
certain polynomials on traces of some other polynomials of a matrix IC^ = S; - for a 

tensor determined by four Stiickelberg helds as 

(6) 

when = ( 1 , 1 , 1 ,— 1 ). A series of arguments presented in [H] (geometrically, we cand consider 
corresponding nonholonomic frame transforms and nonholonomic variables) prove that the parameter 
choice 013 = (a — l)/ 3 ,a 4 = (a^ — a + 1)/12 is useful for avoiding potential ghost instabilities. By 
frame transforms, we can £x 

cx gfiv (7) 

By explicit computations, we can prove that for configurations 0 . de Sitter solutions with 
effective cosmological constant are possible, for instance, for ansatz of PG type, 

ds^ = U‘^{r, t)[dr + e\/ /(r, t)dt]‘^ + d?r‘^{d9‘^ + sin^ Odip"^) — R^(r, t)dt^. (8) 

In above formula, there are used spherical coordinates labelled in the form = r, = 0, = 

ip, y'^ = t), when the function / takes non-negative values and the constant a = a/(a + l) and e = ±1. 
We can consider bimetric conhgurations (determined as soltuons of the system [5]) , [6]) and (I?!)) with 
Stiickelberg helds parameterized in the unitary gauge as (p^ = t and 0- = rn-, 0- = rn-, 0- = rn-, 
where a three dimensional (3-d) unit vector is dehned as n = (n- = sin 6 cos ip, n- = sin^sinip,!!- = 
cos 9). 

First we note that any PG metric of type ([S]) dehnes solutions both in GR and in MGT. For 
instance, we can extract the de Sitter solution, in the absence of matter, and obtain standard cos¬ 
mological equations with FLRW metric, for a perfect huid source 

Tfj.i' [p{t) + p{t)]Uf,u^ + p{t)g (9) 

where = (0,0,0,—R) can be reproduced for the ehective cosmological constant = [R/a. 

Secondly, it is also possible to express metrics of type ([8]) in a familiar cosmological FLRW form (see 
formulas (23), (24) and (27) in [T4]). 

Finally, in this section, we note that oh-diagonal deformations of such solutions can be constructed 
for diherent classes of zero graviton mass and/or massive theories and various f-modihcations, see 
examples in na and m- 
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3 Generating off-diagonal cosmological solutions 

We now make a crucial assumption that our Universe can be described by inhomogeneous cos¬ 
mological metrics are with Killing symmetry on and, in general, can not be diagonalized by 

coordinate transforms. As a matter of principle, we can consider dependencies on all spacetime coor¬ 
dinates but this require more cumbersome computations, see examples in Refs. |12) . Inhomogeneities 
and anisotropies can be very small, but it is important to hnd certain classes of general solutions for 
generic nonlinear systems and impose at the end certain homogeneity and high symmetry conditions 
by selecting corresponding subclasses of generation and integration functions. Up to general classes 
of frame transforms, we can consider the ansatz 

ds"^ = T]i{r,6)gi{r)dr‘^ + r]2{r,6)g2{r)d6‘^ + u‘^{r,6,(p,t){T]3{r,6,t) X 

hs{r, 6)[dip + riiir, 0)dx'^]^ + rnijr, 6, t)h4{r, 6, t)[dt + {wi{x^, t) + Wi{x^))dx^]‘^}. ( 10 ) 

The values rja are called "polarization" functions, where u is the so-called "vertical", v, conformal 
factor. 

For metrics ina, we can consider off-diagonal N-coefficients labelled N^{x^,y‘^), where (for pa- 
rameterizations corresponding to this class of ansatz) 

Nf = riiir, 9) and Nf = Wi{x^, t) -f Wi{x^). 

The data for the "primary" metric are 

°gi{r) = - hi{wif,°g 2 {r) = a‘^r‘^,h 3 = SV^sin^^, ^4 = ^/\fU^ - 

= e^/fU^/h4,W2 = 0,ni = 0, ( 11 ) 

when the coordinate system is such way hxed that the values /, f/, U in ([8]) result in a coefficient gi 
depending only on r. 

Using nonholonomic frame transforms, we can parameterize the energy momentum sources ([9]) 
and effective ([7]) in the form 

T2 =-+ a-'X2) =-!-^( "T + = ( ’"T + "T)(5S, (12) 

for constant values '"T := Mpf{d^f)~^ and "T = , with respect to N-adapted 

frames ([3]). In general, such sources are not diagonal and may depend on all spacetime coordinates. 
Our assumption is that we prescribe a distribution with one Killing symmetry in a moment of time 
and then hnd further evolution with respect to certain classes of nonholonomic frames. We emphasize 
that hxing such N-adapted parameterizations we can decouple the gravitational hied equations in 
MGTs and construct exact solutions in explicit form. 

Let us outline in brief the decoupling property of the gravitational and matter held equations in 
GR and various generalizations/ modihcations studied in details in Refs. |I2|. That anholonomic 
frame deformation method (AFDM) can be applied for decoupling, and constructing solutions of the 
MGT held equations ([5]) with any ehective source parameterized in the form (IT^ . see details in Refs. 

We label the target oh-diagonal metrics as g = {g^ = gigi, ha = Vaha, A(“) with coefficients 
determined by ansatz flTOl) . In these formulas, there is not summation on repeating indices in this 
formula. We shall use brief denotations for partial derivatives: Si'll) = 'il)*,d 2 'ip = 'll)', Os'll) = 'ip^ and 
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^ 4 '^ = i/j*. Computing the N-adapted coefficients of the Ricci tensors, when R\ = = R\,R^]^ 

and R^k are not trivial, we write ([S]) as a system of nonlinear partial differential equations (PDE): 


+ V'" 

ct)*hl 
< + in* 

dkU 


2( ’"T+ "T), 
2h3h4(”^T+ “T), 
0,f3wi - ai = 0, 

rikOJ^ + WkOO*, 


(13) 


(14) 


for 


= In 


\/\h^4 


,7 := (In 


|h4| ' 


a, 


= §-d.<P, (3 

2 /l 3 



(15) 


In above formulas, we consider the system of coordinates and polarization functions are fixed for 
configurations with gi = g 2 = ^ and nonzero values 0* and h*. 

We can extract solutions for the LC-configurations with zero torsion if the coefficients of metrics 
are subjected to additional conditions; 


w* = Si In VI h^l, ei In 


0 , diWj = djWi and n* = 0. 


(16) 


Step by step, the system of nonlinear PDE fll3l] - ffT^ can be integrated in general forms for any u 
constrained by a system of linear first order equations (fT4j) . see details in na. The explicit solutions 
are given by quadratic elements 


ds^ = e'^^^'^\{dx^f + {dxY]+ „r^M d<f + {dk n) dx’^]'^- + A)dxf. 


for any 


4 ^ rnY -|- «'Y') •^1- ^ ^ mY _j_ 

4 = 4, (9,4)" = 9,4", w, + mi = 9,4/4" = 9, A 


(17) 


We can construct exact solutions even such conditions are not satisfied, i.e. the zero torsion conditions 
are not stated or there are given certain sources in non-explicit form. So, the AFDM can be applied 
to generate both off-diagonal metrics and nonholonomically induced torsions. There are various 
physical arguments for what type of generating/ integration functions and sources we have to chose 
in order to construct realistic scenarios for Universe acceleration and observable dark energy/ matter 
effects. 

Coming back to the properties of general solutions, we note that we can generate new classes of 
solutions for arbitrary nontrivial sources, ”*T + “T ^ 0, and generating functions, ^{x^,t) := e'^ 
and Hk = dkn{x'‘). Resulting target metrics are generic off-diagonal and can not be diagonalized via 
coordinate transforms in a finite spacetime region because, in general, the anholonomy coefficients 
fUV for ([3]) are not zero (we can check by explicit computations). The polarization ^-functions from 
fITT)) are 

VI = e^/gi. V2 = e^/h. V3 = <^>V4 ( + “T), g, = ($*) V( (18) 

We conclude that prescribing any generating functions $(r, 6 *,f), n{r,9), u;{r,9,(p,t) and sources 
™T, "T and then computing A{r, 9, t), we can transform any PG (and, similarly, FLRW) metric g = 
[gt, ha,Wi,ni) in MGT and/or GR into new classes of generic off-diagonal exact solutions depending 
on all spacetime coordinates. Such metrics define Einstein manifolds in GR with effective cosmological 
constants determined by ’"T + "T. With respect to N-adapted frames ([3|) the coefficients of metric 
encode contributions from massive gravity, determined by "T, and matter fields, included in '"T. 

We also note that it is possible to provide an "alternative" treatment of flTT)) as exact solutions 
in MGT. In such a case, we have to define and analyze the properties of fiducial Stiickelberg fields 
(/ii and the corresponding bimetric structure resulting in target solutions g = {gi, ha, Nj) : 
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Let us analyze the primary configurations related to 


= {^-= a{T)pa ^-= a{T)pa ^n-,^- = TK ^), 


when the corresponding prime PG-metric g is taken in FLRW form 

ds^ = a^(dp^/(l — Kp^) + p^{d6‘^ + sin^ 6dp^)) — dr'^. 


A corresponding fiducial tensor (jS]) is computed 

o np^ 

Yjuydu^u- = + p^idO"^ + sin^ 9dp^) + 2Hpdpdr — ( „ „ — H‘^p‘^)dT‘^], 

- a^ 

where the coefficients and coordinates are re-defined in the form r ^ p = ar/a{T) and t ^ t = nt, for 
K = 0, ±1; K is an integration constant; H := din a/dr and the local coordinates are parameterized 
in the form x- = p,x- = 9, y- = p,y- = r. 

At the next step, for a target metric g = g^^ and frames e^, = e^da, we write 



9ij + Afj gab 

N^hab ■ 

for e - = 


Nfhab 

hab 


_ 0 e/ _ 


The values 


Pij = = e^5ij = diag[pigi], hah = = diag[paha], Nf = dk n, Nf = dk A 

are related algebraically to data fllSp resulting in off-diagonal solutions fll7p . Then, to work out 
the the "target" Stiickelberg fields we compute 0^ = (/L with being inverse to e^-, and the 

fiducial tensor 

= e 2e/E,a5. 

If the prime value carries information about two constants k and a, a target tensor is 
associated to off-diagonal solutions and encodes data about generating and integration functions 
and via superpositions on possible Killing symmetries, on various integration constants. Similar 
constructions were elaborated for holonomic and nonholonomic configurations in GR, see |16) . 

In the framework of MGT, two cosmological solutions g and g related by nonholonomic defor¬ 
mation^ are characterised respectively by two invariants 


fay ^ ^gaPQjaQ^^y ja/3 _ g“/3e„ 


The tensor does not contain singularities because there are not coordinate singularities on horizon 
for PG metrics. It should be emphasized that the symmetry of is not the same as that of 
and the singular behaviour of 1“^ depends on the class of generating and integration functions we 
chose of constructing a target solution g. 

In GR, MGTs and/or Einstein-Finsler gravity theories na, off-cosmological solutions of type (Ha 
were found to generalized various models of Biachi, Kasner, Godel and other universes. There are 
known locally anisotropic black hole and wormhole, in general, with solitonic background solutions, 
see [min]. For instance, Bianchi type anisotropic cosmological metrics are generated if we impose 
corresponding Lie algebra symmetries on metrics. It was emphasized in [TT] that "any PG-type 
solution in general relativity (with a cosmological constant) is also a solution to massive gravity." Such 

^involving not only frame transforms but also deformation of the linear connection structure when at the end there 
are imposed additional constraints for zero torsion 








a conclusion can be extended to a large class of generic off-diagonal cosmological solutions generated 
by effective cosmological constants but it is not true, for instance, if we consider nonholonomic 
deformations with nonholonomically induced torsion like in metric compatible Finsler theories. 

Finally, we note that the analysis of cosmological perturbations around an off-diagonal cosmo¬ 
logical background is not trivial because the hducial and reference metrics do not respect the same 
symmetries. Nevertheless, fluctuations around de Sitter backgrounds seem to have a decoupling limit 
which implies that one can avoid potential ghost instabilities if the parameter choice is considered 
both for diagonal and off-diagonal cosmological solutions, see details in |18) . This special choice also 
allows us to have a structure ~ at list in N-adapted frames when the massive gravity effects 
can be approximated by effective cosmological constants and exact solutions in MGT which are also 
solutions in GR. 


4 Examples of off-diagonal solutions with solitonic configura¬ 
tions 

We now consider three examples of off-diagonal cosmological solutions with solitonic modihcations 
in MGT and (with alternative interpretation) GR. Two and three dimensional solitonic waves are 
typical nonlinear wave conhgurations which can be used for generating spacetime metrics with Killing, 
or non-Killing, symmetries and can be characterised by additional parametric dependencies and 
solitonic symmetries. Moving solitonic conhgurations can mimic various types of modihed gravity 
dark energy and dark matter effects with nontrivial gravitational vacuum, polarization of constants 
and additional nonlinear diagonal and off-diagonal interactions of the gravitational and matter helds. 


4.1 One soliton solutions 

We shall for simplicity work with a nonlinear radial (solitonic, with left s-label) generating func¬ 
tion 

$ = ^$(r, t) = darctane^'^*'^”'"*^"'"'^® (19) 

and a; = 1, we construct a metric 


ds2 _ de^) + 




4 ( ™T+ “T) 


h3{r,e)d^p^ 


{dt *$)' 


(- mj s^2 


hi{r,t)[dt + {dr A)drY, ( 20 ) 


In this metric, for simplicity, we hxed n{r,6) = 0 and consider that A{r,t) is dehned as a solution 
of ^<|)*/ = dr A and ha are given by PG-data ffTTj) . The generating function flT^ . where 

= (1 — for constants g, go, v, is just a 1-soliton solution of the sine-Gordon equation 

+ sin = 0. 


For any class of small polarizations with rja ~ 1), we can consider that the source ( -|- "T) 

is polarized by when hs ~ hs and ~ h 4 ( with an off-diagonal term dr A 

resulting in a stationary solitonic universe. If we consider that {d^f)~^ = in flT^ . we can 

model /-interactions of type ([T]) via off-diagonal interactions and "gravitational polarizations". 

In absence of matter, = 0, the off-diagonal cosmology is completely determined by "T 
when '^<|) transforms p into an anisotropically polarized/variable mass of solitonic waves. Such 
conhgurations can be modelled if "T. If ^ “T, we generate cosmological models 

determined by distribution oh matter helds when contributions from massive gravity are with small 
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anisotropic polarization. For a class of nonliolonomic constraints on <F and il) (which may be not of 
solitonic type), when solntions (I2U]) are of type (ITU)) with r^Q, ~ 1 and ni,Wi ~ 0, we approximate 
PG-metrics of type ([8]). 

Hence, by an appropriate choice of generating fnnctions and sonrces, we can model eqnivalently 
modified gravity effects, massive gravity contribntions or matter field confignrations in GR and MGT 
interactions. For well defined conditions, snch confignrations can be stndied in the framework of some 
classes of off-diagonal solntions in Einstein gravity with effective cosmological constants. 

4.2 Three dimensional solitonic anistoropic waves 

More sophisticate nonlinear gravitational and matter filed interactions can be modeled both in 
MGTs and GR if we consider more general classes of solntions of effective Einstein eqnations. 

For instance, off-diagonal solitonic metrics with one Killing symmetry can be generated, for 
instance, if we take instead of fll9p a generating fnnctions ^$(r, 6*,f) which is a solntion of the 
Kadomtsev-Petivashvili, KdP, eqnations |19| . 

when solntions indnce certain anisotropy on In the dispersionless limit e —)■ 0, we can consider 
that the solntions are independent on 6 and determined by Bnrgers’ eqnation 

^<1* = 0. 

Snch solntions can be parameterized and treated similarly to fl2Up bnt with, in general, a nontrivial 
term [de A)d9 after h^, when = A* and *<!>'/ = A'. Similar metrics were constrncted 

for Dirac spinor waves and solitons in anisotropic Tanb-NUT spaces and in five dimensional brane 
gravity which can be encoded and classihed by corresponding solitonic hierarchies and geometric 
invariants, see czi- Here we proved that AFDM can extended to generate inhomogeneons off- 
diagonal cosmological solitonic solntions in varios MGTs. 


4.3 Solitonic waves for a nontrivial vertical conformal t;—factor 


The cosmological solntions we look are also with three dimensional solitons bnt for a n-factor 
as in oa, For instance, we consider solitons of KdP type, when u = u}{r,(p,t)j when = r,x‘^ = 
9, = ip,y‘^ = t, for 

± + {dtu + uu* + eo;”*)* = 0, (21) 

In the dispersionless limit e ^ 0, we can consider that the solntions are independent on angle ip and 
determined by Bnrgers’ eqnation 

(h* +(1; (!;• = 0. 

The conditions fll4p impose an additional constraint 

eica = cj* + wi{r, 9, (p)uj* + ni(r, 9)u'^ = 0. 


In the system of coordinates when u' = 0, we can fix ^2 = 0 and n 2 = 0. For any arbitrary generating 
fnnction with LG-confignration, ^{r,9,t), we constrnct exact solntions 


ds2 _ d9^) + 




4 I" mY oX) 


h^{r,9)dyA 




I T + “T) $2 


hi{r,t)[dt+{dr A)dr]‘^, (22) 


^In a similar form, we can construct various types of vacuum gravitational 2-d and 3-d configurations characterized 
by solitonic hierarchies and related bi-Hamilton structures, for instance, of KdP equations with possible mixtures with 
solutions for 2-d and 3-d sine-Gordon equations etc, see details in Ref. El¬ 


io 





which are generic off-diagonal and depend on all spacetime coordinates. Such stationary cosmological 
solutions are with polarizations on two angles 6 and ip. Nevertheless, the character of anisotropies is 
different for metrics of type (1201) and (l22|) . In the third class of metrics, we obtain a Killing symmetry 
on only in the limit a) —)■ 1, but in the first two ones, such a symmetry exists generically. For 
fl22|) . the value $ is not obligatory a solitonic one which can be used for additional off-diagonal 
modifications of solutions and various types of polarizations. 

We can provide a physical interpretation of [22] which is similar to [20] if the generating and 
integration functions are chosen to satisfy the conditions ~ 1 and ni,Wi ~ 0, we approximate 
PG-metrics of type ([8]). 

In a particular case, we can use a conformal n-factor which is a l-solitonic one, i.e. 

Cj oj{r,t) = 4 arctan , 

where cx^ = (1 — and constants q, Qq, v, defines a 1-soliton solution of the sine-Gordon equation 

u — u + sinca = 0. 

Such a soliton propagates in time along the radial coordinate. 

We conclude that solitonic waves may mimic both particle type configurations as dark matter 
and encode certain hidden dark energy and off-diagonal gavitational and matter field interactions. 

5 Reconstruction mechanism for off—cosmological solutions 

We now consider a reconstruction mechanism with distinguished off-diagonal cosmological effects 
n by generalizing some methods elaborated for f{R) gravity in [3j. The main idea was to present 
the MGT actions (with zero or non-zero gravitational mass) as sums of actions in GR and certain 
effective ideal fluid contributions with parameters defined by nontrivial /i and /-deformations. The 
reconstruction method was developed for such formulations which lead to cosmology with cyclic 
evolution. Then, it was proven that the ekpyrotic scenario may be also realized for MGTs and that 
it is possible to reconstruct models of f{R) gravity which induces little rip cosmology. 

In our approach, we work with cosmological generic off-diagonal metrics and generalized connec¬ 
tions. We can always chose such generating functions and parameters of solutions that off-diagonal 
contributions are small, and the torsion is constrained to be zero, beginning a certain fixed moment 
of time. Nevertheless, we can not neglect for such classes of solutions possible nonlinear effects deter¬ 
mined by generating functions and effective sources resulting in diagonalized modifications. MGTs 
were elaborated as realistic alternatives for unified description of inflation with dark energy when 
cosmological scenarios encode information on various massive and zero-mass gravitational modes. A 
crucial question is if and how such constructions have to be elaborated for generic off-diagonal cos¬ 
mological spaces. This is not only a geometric problem for generalizing the reconstruction formalism 
for inhomogeneous and locally anisotropic cosmological theories. It is connected to a very important 
question on equivalent modeling of cosmological scenarios for different MGTs in the framework of 
GR with nonholonomic and off-diagonal nonlinear interactions. 

Any cosmological solution in massive, MGT and/or GR parameterized in a form flTU]) (in partic¬ 
ular, as (12n|) and (I22|) I can be encoded into an effective functional (0]) when 

This allows us to work as in MGT; the conditions daf{R) = 0 ii R = const simplify substantially 
the computations. Nevertheless, the contributions of the parameter p and effective potential U are 
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included as functional dependencies of / and i?.For D ^ V, we get nonlinear modifications both in 
diagonal and off-diagonal terms of cosmological metrics. 

The starting point of our approach is to consider a prime flat FLRW like metric 

ds"^ = a^{t)[{dx^y + {dx^Y + — dt'^, 


where t is the cosmological time. In order to extract a monotonically expanding and periodic cosmo¬ 
logical scenario, we parameterize In |a(t)| = ifo^+d(t) for a periodic function a(t-|-r) = cos(27rt/r), 
where 0 < < Hq. Our goal is to prove that such a behavior is encoded into off-diagonal solutions 

of type (I 20 |) - fl 22 |) . 

We write FLRW like equations with respect to N-adapted (moving) frames ([3]) for a generalized 
Hubble function H, 

3H^ = Sirp and 3H^ + 2e^H = —Sirp. 

Using variables with consider a function H{t) when e^H = dfH = H*. 

It should be noted that the approximation 64 —)■ dt is considered "at the end" (after a class of 
off-diagonal solutions was found for a necessary type connection, which allowed to decouple the 
equations) and the generating functions and effective sources are constrained to depend only on t. 
The energy-density and pressure of an effective perfect fluid are computed 

p = (8;r)-'[(5„/)-‘(i/(fl) + 3ffe4(SK/))-3e4i/] (23) 

= (87r)-'13gIti Tl/l - 3H-] = (87r)-‘|% In ^J\f^U + f(R)\ - 3H-], 

P = -(8ir)-'[(5K/)“'(i/(fl) + 2i/e4(5R/)+e 4 e 4 (SK/))+ 64 //] 

= (87r)-'[SjjIn Tl/l + H-] = -(87r)-‘19, In + + H']. 

We emphasize that, in general, such values are defined with respect to nonholonomic frames. The 
effect of nonlinear deformations encoding physical data f, R) is preserved also in diagonal 

contributions for 64 —)■ 5^. 

The equation of state, EoS, parameter for the effective dark fluid encoding MGTs parameters is 
defined by 

^ p f + 2H‘agf fM +/(g) + 2g»ag/ 

P f- 6H>dff f M + f(R) - 6H-dgf' 

when the corresponding EoS is 

p = —p — 

and W(t) is computed, for simplicity, for a configuration of "target" Stiickelberg fields (p^' = 
when a found solution is finally modelled by generating functions with dependencies on t. The 
components of 0^ are computed with respect to N-adapted frames. 

Taking a generating Hubble parameter H{t) = Hot -|- Hi sin ut, for cu = 27r/r, we can recover the 
modified action for oscillations of off-diagonal (massive) universe (see similar details in [21), 


f{R{t)) = 6uHi 



4 cosut{Ho + Hi sin cut)] exp[iLot -|- 


u 


sin ut]. 


(25) 


Here we note that we can not invert analytically to find in explicit form R, or any nonholonomic 
deformation to R with respect to general N-adapted frames. Nevertheless, we can prescribe any 
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values of constants Hq and Hi and of oo and compute effective dark energy and dark matter oscillating 
cosmology effects for any off-diagonal solution in massive gravity and/or effective MGT, GR. To 
extract contributions of fi we can fix, for instance, f{R) = R = R and using ([1]) and ([2]) we can relate 
f{R{t)) and respective constants to certain observable data in cosmology. 

The MGT theories studied in this work encode, for respective nonholonomic constraints, the 
ekpyrotic scenario which can be modelled similarly to f{R) gravity. A scalar field is introduced 
into usual ekpyrotic models in order to reproduce a cyclic universe and such a property exists if 
we consider off-diagonal solutions with massive gravity terms and/or /-modifications. The main 
idea is to develop the reconstruction techniques for the scalar-tensor theory using the AFDM, with 
nonholonomic off-diagonal metric and linear connection deformations. Working with general classes 
of off-diagonal cosmological solutions, the problem is to state the conditions for generating and 
integration functions when a corresponding ekpyrotic scenario will "survive" for certain constraints 
and in respective limits. 

Let us consider a prime configuration with energy-density for pressureless matter pm, for radiation 
and anisotropies we take respectively pr and p^ for radiation and anisotropies, k, is the spatial 
curvature of the universe and a target effective energy-density p fl23|) . A FLRW model can be 
described by 


3ff2 = 8x[^ + ^ + ^-4 + p]. 

a-^ H 

We generate an off-diagonal/massive gravity cosmological cyclic scenario containing a contracting 
phase by solving the initial problems if tc > 1, see fl2T|) . A homogeneous and isotropic spatially flat 
universe is obtained when the scale factor tends to zero and the effective /-terms (massive gravity 
and off-diagonal contributions) dominate over the rests. In such cases, the results are similar to 
those in the inflationary scenario. For recovering fl25p . the ekpyrotic scenario takes place and mimic 
the observable universe for f ~ 7r/2ci; in the effective EoS parameter 


tc ~ —1 + sinwf/ScniLi cos^ ut S> 1. 


This allows us to conclude that in massive gravity and/or using off-diagonal interactions in GR cyclic 
universes can be reconstructed in such forms that the initial, flatness and/or horizon problems can 
be solved. We can compute possible locally anisotropic and inhomogeneous small contributions for 
self-consistent models with nonholonomic frames. 

In the diversity of off-diagonal cosmological solutions which can constructed using above presented 
methods, there are cyclic ones with singularities of the type of big bang/ crunch behaviour. This 
is still a largely unexplored area both for the geometric methods of constructing exact solutions of 
PDEs and recovering procedures for certain "preferred" fundamental physical objects compatible 
with experimental data. Choosing necessary types generating and integration functions, we can 
avoid singularities and elaborate models with smooth transition. Using the possibility to generate 
nonholonomically constrained /-models with equivalence to certain classes of solutions in massive 
gravity and/or off-diagonal configurations in GR, we can study in this context, following methods 
in mm, big and/or little rip cosmology models, when the phantom energy-density is modelled by 
off-diagonal interactions. We note that such nonsingular models for dark energy were proposed as 
alternatives to ACDM cosmologies, [20] ■ Using corresponding classes of generating functions, we can 
reproduce the sceanaros with a phantom scalar modeling a little rip in the framework of AFDM and 
nonholonomic MGTs. We omit such considerations in this article (see a summary of such construction 
and further developments in jl5|i. 
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6 Conclusions 


In this paper we deal with new classes cosmological off-diagonal solutions in massive gravity 
with flat, open and closed spatial geometries. These solutions can be systematically constructed 
for various types of modified gravity theories, MGTs, and in general relativity, GR. We applied an 
advanced geometric techniques for decoupling the field equations and constructing exact solutions 
in massive and zero mass f{R) gravity, theories with nontrivial torsion and noholonomic constraints 
to GR and possible extensions on (co) tangent Lorentz bundles. The so-called anholonomic frame 
deformation method, AFDM, was elaborated during last 15 years in our works Ha E], where a 
number of examples of off-diagonal solutions and new applications in gravity and modern cosmology 
were considered. 

A very important property of such generalized classes of cosmological solutions is that they 
depend, in general, on all spacetime coordinates via generating and integration functions and con¬ 
stants. They describe certain models of inhomogeneous and locally anisotropic cosmology with less 
clear physical meaning and possible physical implications [5]. After some classes of solutions were 
constructed in a most general form, we can impose at the end additional nonholonomic constraints, 
cosmological approximations, extract configurations with a prescribed spacetime symmetry and/or 
dependence on certain mass parameters, consider asymptotic conditions etc. Thus, our solutions 
can be used for elaborating homogeneous and isotropic cosmological models with arbitrary spatial 
curvature, to study generalized Killing and non-Killing, with possible nonholonomically deformed 
(super) symmetries [TTl [T6] and to study "non-spherical" collapse models of the formation of cosmic 
structure such as stars and galaxies (see also [S]). 

Another aspect of the AFDM is that if we work only with cosmological scenarios for diagonalizable 
metrics, there are possibilities to discriminate the massive gravity theory form the /-gravity and/or 
GR. For diagonal metrics depending on a time like coordinate, we can formulate mathematical 
cosmology problems for certain nonlinear systems of ordinary differential equations with general 
solutions depending on integration constants. Identifying such a constant, for instance, with a 
graviton mass parameter, we do not have much possibilities to mimic a number of similar effects in 
GR for different MGTs. Following only such a "diagonal" approach, we positively have to modify 
the GR theory in order to explain observational data in modern cosmology and elaborate realistic 
quantum models of massive gravity. 

We now discuss a new and important feature of the off-diagonal anisotropic configurations which 
allows us to model cosmic accelerations and massive gravity and/or dark energy and dark matter 
effects as certain effective Einstein spaces. Having integrated such system of nonlinear partial differ¬ 
ential equations, PDEs, for a large class of such solutions, we can put and motivate such questions; 
May be we do not need to modify radically the GR theory but only to extend the constructions to 
off-diagonal solutions and nonholonomic systems and try to apply this in modern cosmology? Gould 
we explain observational data in modern cosmology via nonlinear diagonal and/or off-diagonal in¬ 
teractions with non-minimal coupling for matter and/or different phases of massive and zero mass 
gravity. This is a quite complicated theoretical and experimental problem and the main goal of this 
and our recent papers cited in Refs. da [15] was to analyze such constructions from the viewpoint 
of massive gravity theory when off-diagonal effects can be alternatively explained to other types of 
gravity theories. 

The reconstruction procedure for cosmological models with non-minimally coupled scalar fields 
evolving on a flat FLRW background and in different MGTs was studied in HE]. In this work, 
we elaborated a reconstruction method for the massive gravity theory which admits an effective 
off-diagonal interpretation in GR and /-modified gravity with cyclic and ekpyrotic universe solu¬ 
tion. We concluded that the expansion can be around the GR action even if we admit a nontrivial 
effective torsion. For zero torsion constraints, it is possible to construct off-diagonal cosmological 
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models keeping the approach in the framework of the GR theory. We further investigated how our 
results indicate that theories with massive gravitons, with possible /-modihed terms and off-diagonal 
interactions may lead to more complicated scenarios of cyclic universes. Following such nonlinear 
(off-diagonal) approaches, the ekpyrotic (little rip) scenario can be realized with no need to introduce 
additional fields (or modifying gravity) but only in terms of massive gravity or GR. Another inter¬ 
esting constructions can be related to reconstruction scenarios of f{R) and massive gravity theories 
leading to little rip universes both in locally anisotropic and isotropic variants. Finally, we note that 
the dark energy for little rip models present an example of non-singular phantom cosmology. 
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